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Abstract 

We give a physical interpretation of the recently demonstrated nonconservative nature of interatomic forces in current-carrying 
nanostructures. We start from the analytical expression for the curl of these forces, and evaluate it for a point defect in a current- 
carrying system. We obtain a general definition of the capacity of electrical current flow to exert a nonconservative force, and thus 
do net work around closed paths, by a formal noninvasive test procedure. Second, we show that the gain in atomic kinetic energy 
over time, generated by nonconservative current-induced forces, is equivalent to the uncompensated stimulated emission of direc- 
tional phonons. This connection with electron-phonon interactions quantifies explicitly the intuitive notion that nonconservative 
forces work by angular momentum transfer. 



Introduction 

Electron-nuclear interactions lie at the heart of the transport 
properties of nanoscale conductors. Even in the limit of elastic 
(phonon-free) conduction, the nature and positions of nuclei in a 
nanojunction determine the external potential, experienced by 
current-carrying electrons, and, together with electron-electron 
interactions, determine the current-voltage spectrum of the 
system. Allowing nuclei to respond to current-induced forces 



introduces two additional elements: Current-driven displace- 
ments and Joule heating. Current-induced forces arise, funda- 
mentally, through momentum transfer from the electron flow to 
nuclei, and are familiar from the field of electromigration [1]. 
An alternative, but fundamentally equivalent, way to think 
about them is as nonequilibrium corrections to interatomic 
bonding forces. In considering these forces, it is often conveni- 
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ent to adopt the Bom-Oppenheimer approximation, adapted to 
the nonequilibrium conditions in a nanoconductor: We think of 
nuclei as heavy, slow classical particles, and imagine that, as 
nuclei move, electrons always remain in the steady state appro- 
priate for the given set of instantaneous nuclear positions; we 
then calculate the force on a nucleus exerted by the mean elec- 
tron density in the system (including appropriate Pulay correc- 
tions, if an incomplete electronic basis is used) [2-5]. Joule 
heating, on the other hand, is due to the finite mass of nuclei, 
and results from the recoil of nuclei in inelastic collisions with 
electrons [6-10]. The combined effect of the two is the driving 
force behind electromigration-type phenomena [2,6,7]: Current- 
induced forces modify atomic migration barriers; together with 
local heating, this results in thermally activated current-induced 
atomic rearrangements, or even failure. 

Recently, a new and rather different aspect of current-induced 
forces has received attention: Their nonconservative character, 
and their resultant ability to do net work on individual atoms, or 
groups of atoms, around closed paths [1,11-14]. The practical 
consequences of this mechanism for sustained energy transfer 
from electrical current into atomic motion are only just starting 
to be explored. Two aspects of the effect that are of immediate 
interest are its capacity to drive an atomic-scale motor, and its 
possible potent role as a cause for dramatic mechanical failure 
[12-14]. Indeed, the notion that current can drive rotary motion, 
under appropriate conditions, is highly intuitive [15,16], and is 
increasingly being seen as a common, rather than rare, effect in 
nanoscalc conductors [17]. The essential physics behind the 
nonconservative component of current-induced forces is that of 
a waterwheel driven by a flow [12]. For quantum-mechanical 
electron flow interacting with classical nuclei, the effect is 
quantified precisely by an analytical result for the curl of these 
forces [12,18]. Yet, there is a gap to be bridged between the 
formal result and the intuitive physics. 

The aim of the present short paper is to bridge this gap, and 
extract explicitly the gas-flow picture of nonconservative 
current-induced forces. We then make a second coimection, by 
showing that the work done by these forces around closed paths 
is equivalent to the uncompensated stimulated emission of 
directional phonons, characterised by the sign of their angular 
momentum. This second result will close the gap between the 
nonconservative effect and the more familiar fundamental 
physics of electron-phonon interactions. 

Results and Discussion 
The gas-flow picture 

Under steady-state conditions, in the absence of phonons, the 
electronic properties of a nanoscale conductor are parametric 
functions of the classical nuclear positions. So too are the 



current-induced forces on the nuclei. The nonconservative 
component of these forces is characterised by the generalised 
curl expression [12,18] 



BFjiQ) 8Fj(Q) 
8Qi dQj 



: 4TzbnTr{FjiQ)Diii,Q)FjiQ)Ap{Q)}. (i) 



Here, and throughout the paper, we work in the small-bias limit. 
For our present purposes, we work with noninteracting elec- 
trons and zero magnetic fields. 

In Equation 1 , Q denotes a collection of generalised classical 
coordinates {Qj}; F^Q) is the current- induced force on co- 
ordinate Qi, D{\i,Q) = 5\\i - Hg(Qy\ is the operator for the elec- 
tronic density of states, evaluated at the Fermi level |i; Ap(0 is 
the nonequilibrium part of the one-electron density matrix (that 
is, the difference between the steady-state current- carrying 
density matrix and the equilibrium density matrix); finally, 
Fj{Q) is the force operator 



FiiQ)- 



(2) 



where Hp(Q) is the one-electron Hamiltonian, as a parametric 
function of the classical degrees of freedom Q. 

Equation 1 is intriguing but lacks transparency. We now probe 
its physical content as follows. We immerse a point defect, with 
position vector R = {X,Y,Z) and a scattering potential v(R) = 
CJ(f-R), in the electron flow. Here, r = {x,y,z) denotes the 
electron position. We shall use the defect to directly measure 
the ability of the flow to exert a nonconservative force. 



We have 



F(R) = -CVR5(f - R) = C Vr8(f - R). 



(3) 



Below, we deliberately treat C as small and work to the lowest 

nontrivial order in this parameter. Taking the trace in 
Equation 1 in the continuum r-representation. 



(Vr xF(R))2 = 47ilmj drdr'F^(r,R)£)(r,r',R,n) 

xF},(r',R)Ap(r',r,R). 



(4) 



We also have 



m \ ov 



(5) 
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(VrXj(r,R)), =— Im 
m 



'^5^Ap(r,r',R)^ 



dydx' 



(6) 



/r=r' 



where v = x,y,z and j(r,R) is the electron particle-current 
density. 

Then, to lowest order in C, we obtain 



Finally, recognising 



(10) 



as the electron-wind force in electromigration, we have 

VrxF(R) = (R). (11) 



VrXF(R): 



2%mC 



(VrX[D(r,^)j(r)]}^^j^, (7) 



where D(r,\i.) = D(r,r,|a,) and j(r) are the local density of states 
at the Fermi level and the particle-current density in the absence 
of the point scatterer. Therefore, Equation 7 defines an intrinsic 
property of the current flow: Its "curl-generating" capacity. The 
point defect above serves as a noninvasive test particle that 
probes this property. One can picture situations in which this 
intrinsic curl vanishes but yet there still is a nonzero curl to 
higher order in the coupling between the scatterer and the elec- 
trons. One example of this interesting possibility would be an 
atom weakly bonded to a structure, where the weak bonding 
enables current to flow through the atom in the first place. 

Equation 7 is completely general and makes no assumptions 
about the nature and structure of the conductor. We will now 
simplify it further as follows. We assume that all electronic 
properties vary slowly in space, and we attribute a local Fermi 
momentum, /»(r), to the electrons. We now have a semiclassical 
gas flow in a, locally, jellium-like environment. Next, we 
observe that for jellium 



h 



p{r)/m 



(8) 



r=R 



is the scattering cross section of the defect. In the assumed, 
locally free-electron-like medium, in 3D, to lowest order in C, a 
is a constant (since D{r,[i) x p{r)), and is equal also to the 
transport cross section, Otr- Hence, 



Therefore, we have shown from first principles that the point of 
departure. Equation 1, is an algebraic statement of Sorbello's 
thought experiment [1] to prove that this force is, in general, a 
nonconservative force. We have shown, further, that the key 
quantity responsible for this property is the curl of the local 
electron momentum-current density. 

Nonconservative work as directional plionon 
emission 

We will now relate the nonconservative current-induced forces 
on atoms to the intuitive idea of the waterwheel. To this end, we 
will show that the work done by these forces around closed 
paths corresponds exactly to the stimulated generation of direc- 
tional phonons, characterised by their angular momentum. 
Consider a flux of electrons travelling through an elastic 
medium. Intuitively, we expect the flux to preferentially emit 
phonons with momentum parallel to the electron flow. This 
becomes evident, when we consider the setup in Figure 1 . 




Vr X F(R) = a^, { x [;.(r)j(r)]}^^j^ . (9) 

Thus, (the curl of) the force on the test particle is proportional 
to (the curl of) the local momentum flux of the flow, [p(r)j(r)], 
with a constant of proportionality Ot^. At this stage, quantum 
mechanics has all but disappeared from the problem: We have a 
classical interaction between a, generally, spatially nonuniform 
steady gas flow and an elastic scatterer in its path. 



V 

Figure 1 : An electrode-junction-electrode system In the Landauer 
picture. The details are discussed in the text. 

A nanostructure is connected to two electrodes, each in turn 
connected to its own battery terminal serving as a particle reser- 
voir. The left reservoir injects right-travelling electrons with 
electrochemical potential 
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Ill = ii + eV / 2, 



with all other commutators equal to zero. The inverse relations 
are 



and the right reservoir injects left-travelling electrons with elec- 
trochemical potential 



\1r = \i-eV 12, 



(13) 



X 



2VM(o V + 



(18) 



where V is the applied bias. 

The current comes from the energy window between \ir and H£, 
where there are partially-populated electron states. A right- 
travelling electron in this energy window can emit a right- 
travelling phonon, and scatter into a left-travelling state. But the 
reverse process is suppressed, due to the population imbalance 
between the two sets of electron states. Hence, we expect a 
directional preference of the emitted phonons. Since the rate of 
stimulated emission increases with growing phonon population, 
we might expect the process to have the ability to feed on itself. 
We will show below that this is the physical origin of the 
nonconservative current-induced forces. 

Consider two independent generalised oscillator coordinates X 
and Y, with the same angular frequency, co. Here, X and 7 could 
be two individual atomic degrees of freedom (not necessarily of 
the same atom), or they could be collective normal modes. In all 
cases, we assume, as our starting point, that the modes X and Y 
describe standing waves. The Hamiltonian for the two oscil- 
lator degrees of freedom is 



(14) 



where Px and Pj are the corresponding canonical momenta, and 
M is a mass-like parameter. 

We now form new modes, which we label by {+) and (-), with 
annihilation operators 



1 Mfi) 



(x-iY)- 



i 1 



2 V Mho) 



(Px-A) (15) 



A ^1 
2 

They obey 



i^x+iiy 



1 



2 VM^to' 



(16) 



[i^,4] = l = [i,/f], 



(17) 



Y =- 



1 h 



2i VMoo 



(19) 



Pj^ =^Vm^co {a++A_ -ij (20) 

Py = ^ VmAco (i+ +^-A_- Jfj . (21) 
The new modes diagonalise the Hamiltonian as follows 



^osc 



iV, -l---l-.^_-l-- 
^ 2 2 



(22) 



where 



(23) 



The physical significance of the new modes can be seen by 
considering the angular momentum 



Lz=XPY-YPj,=h(N^-N_). 



(24) 



Thus, the "anticlockwise" mode (-I-) carries positive angular 
momentum (in the direction, perpendicular to the abstract X-Y 
plane) and the "clockwise" mode (-) carries negative angular 
momentum. By couphng these two directional phonon modes to 
electrons, we will see that the electron current pumps energy 
into one, while damping the other. 

The coupling between electrons and phonons is described by 
scattering theory. The unperturbed, phonon-free state of the 
current-carrying electrons is that of the usual Landauer picture 
[5]. In this picture, electrons in the phonon-free steady state are 
described by two sets of stationary one-electron 
Lippmann-Schwinger scattering states. One set, with 
energies {£/}, originate from the left electrode, and are 
scattered elastically at the junction, with partial backscattering 
into the left electrode and partial transmission into the right 
electrode; the other set, with energies {Er}, originate 
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from the right electrode, with partial backscattering into the 
right electrode and partial transmission to the left. In the 
absence of bound states (which we assume here), the two form a 
complete orthonormal set. The left reservoir populates states /, 
in a grand-canonical ensemble with electrochemical potential |X£ 
(and a chosen electronic temperature); the right reservoir popu- 
lates states r with electrochemical potential It is convenient 
to defme the density of states operators [5] 



with 



I 



(25) 



where /; =fiiEi) and = fR(Er) are the Fermi-Dirac distribu- 
tions for electrons originating from the two respective reser- 
voirs. 

We now deliberately suppress the spontaneous phonon emis- 
sion. Formally, we work in the classical limit and set {N± +\)~ 
N±. Then, counting all possibilities for a and p above and 
observing the selection rules, setting the electronic temperature 
to zero, and ignoring variations in the electronic properties over 
energies in the region of hm or eV, we get 



N+ nh 
2 M 

Mco 



(30) 



DLiE) + Dj,iE)=DiE). 



(26) 



where 



Here, we include spin in Di(E) and Dj^{E). From Equation 12, 
Equation 13 and Equation 25, the nonequilibrium part of the 
one-electron density matrix, in the linear bias regime, is given 
by 



Ap = ^[4(M)-^«(p)]. 



(27) 



The system of two oscillators, we describe by an unperturbed 

density matrix that commutes with A'+. Electrons and oscil- 
lators are then coupled by an interaction of the generic form 



(28) 



where c„ and c\ are fermionic annihilation and creation oper- 
ators. Then, for the resultant rates of change of the occupancies, 
7V±, of the anticlockwise and clockwise modes, to lowest order 
in V, we obtain 



2M(o 

xZ Z /a(l-/p)KVpl(^x + i^y)IV„>l 

a=l,r P=/,r 

X 5(£'p +fi(0) 



-N+ 



(29) 



2Mfi) 

xZ Z /a(l-/p)KVpl(^X±i^F)l¥a>l 
a=/,rp=/,r 



Tj, =Ti:{DL(y.)Fj,Dj^i^)Fj,} 

+ nr{Dj^(li)Fj^Dj^{lx)Fj^} (31) 
+ Tr{Dji{^)Fj^Djii^)Fj^} 

ry=Tr{4(^)/v4(^i)^y} 

+ 2Tr{4(H)/VD^(M)Fj.} (32) 

+ TT{DR(ii)FYDKiii)FY}. 

This is our final result. Equation 30 displays precisely the 
picture from [12]. Mode (+) experiences a damped driven 
motion. The damping (first term) is due to the ordinary elec- 
tronic friction experienced by the two independent modes and 
y(each of which carries half of the energy of mode (+)). This 
Iriction is due to phonon absorption by electrons, and is present 
even at zero current. The driving term (the second term) comes 
solely from the current. It anti-drives mode (-). 

To extract explicitly the curl of the effective driving force 
acting on the 2D oscillator, for the rate at which this force does 
work on mode (+) we write 



=(VxF)2 7l4 



2 to 
271 ' 



(33) 



where (7+ = N+ hm = May^A^ is the energy of the mode and 
is the mean square displacement in the X-F plane. Comparison 
with the second term in Equation 30 gives 

(V X F)z = 4neV]mTT{Dji (34) 
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In view of Equation 26 and Equation 27, this is equivalent to 
Equation 1 [13]. The difference is that in [12,18] this result was 
derived in the mixed quantum-classical framework of Ehrenfest 
dynamics. Herein, it has been obtained formally exactly from 
quantum-mechanical electron-phonon coupling: The work done 
by nonconservative forces around a closed path is equivalent to 
the stimulated emission of directional, angular-momentum- 
carrying phonons. The opposite sign of the effect for modes (+) 
and (-) originates from the directionality, introduced by the 
current, that breaks the symmetry between clockwise and anti- 
clockwise atomic motion in the X-F plane [17]. The averaging 
over classical trajectories, implicit in the construction of the 
unperturbed phonon density matrix in the present calculation, 
eliminates certain additional forces that become apparent, for 
example, in the treatments of references [13,14]. These add- 
itional forces and their effects present an interesting avenue for 
further work [19]. 

Conclusion 

We have taken the analytical result for the curl of current- 
induced forces, derived in references [12,18], and we have 
related it to two physically transparent ideas. One is the elec- 
tron-wind force on a test particle, determined by the local elec- 
tron momentum flux. The second is the uncompensated stimu- 
lated emission of travelling phonons. The fact that the noncon- 
servative effect is related to stimulated (as opposed to spontan- 
eous) emission explains the remarkable, and practically very 
useful, earlier finding that the nonconservative dynamics of 
atoms under current can be captured already at the level of 
Ehrenfest dynamics. (Ehrenfest dynamics suppresses spontan- 
eous transitions but retains stimulated transitions.) However, at 
that level, the underlying physics of the effect remains some- 
what obscured by the mixed framework. Starting from the 
present internally consistent picture, we see explicitly that this 
novel and interesting effect is like a current-driven waterwheel, 
which works by angular-momentum transfer from the electron 
flow (albeit possibly in an abstract sense, depending on what X 
and y denote). 

Superficially, Equation 30 resembles ordinary Joule heating. 
Indeed, it does constitute a form of directional heating. But 
there is a key difference. Standing, bound phonon modes can 
equilibrate with the current-carrying electrons, at an effective 
phonon population set by the bias [20]. Equation 30 shows 
qualitatively different behaviour. If the second term is positive, 
then, once the bias V is large enough for term 2 to dominate 
term 1, the Equation predicts an exponential growth of the 
energy of the given travelling mode, without equilibration. This 
is the waterwheel effect of reference [12]. Of course, in reality 
this increase cannot continue indefinitely, and the possible 
outcomes form the subject of ongoing research. Pertinent ques- 



tions are concerned with the effects of anharmonicity [12], the 
possible eventual failure of the device [13], and the possible 
current-induced corrections, under appropriate conditions, to the 

conservative part of the harmonic potential (which could lift the 
degeneracy of the modes Xand 7 that form the "waterwheel"). 

It is hoped that nonconservative forces, and the underlying 
mechanism of uncompensated directional phonon generation 
under current, will be useful in tackling not only problems in 
nanotechnology but also in neighbouring areas, such as the 
behaviour of bulk metals under large current densities. 
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